Ground states of dipolar gases in quasi-lD ring traps 
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We compute the ground state of dipoles in a quasi-one-dimensional ring trap using few-body 
techniques combined with analytic arguments. The effective interaction between two dipoles depends 
on their center-of-mass coordinate and can be tuned by varying the angle between dipoles and the 
plane of the ring. For weak enough interactions, the state resembles a weakly interacting Fermi 
gas or an (inhomogeneous) Lieb-Liniger gas. A mapping between the Lieb-Liniger and the dipolar- 
gas parameters in and beyond the Born approximation is established, and we discuss the effect of 
inhomogeneities based on a local-density approximation. For strongly repulsive interactions, the 
system exhibits crystal-like localization of the particles. Their inhomogeneous distribution may 
be understood in terms of a simple few-body model as well as a local-density approximation. In 
the case of partially attractive interactions, clustered states form for strong enough coupling, and 
the dependence of the state on particle number and orientation angle of the dipoles is discussed 
analytically. 

PACS numbers: 67.85.-d, 05.30.Jp, 05.30.Fk 



I. INTRODUCTION 

The creation and study of so-called dipolar quantum 
eases have recently become a major research focus P, 
2]. These gases, such as ultracold atoms with magnetic 
dipole moments (e.g., Crjp or Dy [H)or polar molecules 
(such as KRb @, LiCs @ or RbCs 0), are dominated 
by dipole-dipole rather than van-der-Waals interactions. 
Combining a high level of control with long-ranged and 
strongly anisotropic interactions, dipolar gases offer both 
the simulation of elusive quantum states in condensed- 
matter or nuclear physics and the design of exotic novel 
quantum phases. 

The question naturally arises how a dipolar gas be- 
haves when confined to lower dimensions. This is im- 
portant experimentally, because the lower dimensionality 
may help alleviate the collisional instability toward head- 
to-tail alignment of the dipoles [8, 9], and conceptually, 
given that the enhanced quantum fluctuations in lower 
dimensions give rise to intriguing physics. 

In particular, in strictly one dimension (ID), dipolar 
gases have been shown to exhibit Luttinger-liquid behav- 
ior (see, e.g., [ld]| - fl3l | ) . A finite, quasi-lD transverse con- 
finement may both remove the short-distance divergence 
of the dipolar interactions [l4| and alter the short-range 
s-wave interactions, which may lead to intriguing physics 
like a roton instability toward a density wave. In spite 
of displaying long-range interactions, this system still is 
an effectively homogeneous ID system, so long as the 
dipoles have a common orientation. In a linear geometry, 
the anisotropy enters only once the system ceases to be 
one-dimensional, in which case "zig-zag" chains or other 
higher-dimensional configurations may show up [l5j . 




A dipolar system exhibiting truly ID physics as well as 
anisotropic interactions can be achieved in curved lower- 
dimensional geometries. Nontrivial geometries are in- 
teresting not only in the context of dipolar interactions 
-20|, but also for studying, e.g., persistent currents 
[22|. In the simplest case of a quasi- ID ring trap, 
which can be realized in the context of cold atoms |23l — 
l26l |. the effective ID dipole interaction between two par- 
ticles becomes inhomogeneous in the sense that it ac- 
quires a dependence on their center-of-mass coordinate. 
In this paper, we show how by varying both the over- 
all strength as well as the degree of inhomogeneity, in- 
teresting regimes are found in such a system - such as 
a ID Bose-gas-like phase with center-of-mass dependent 
interaction 7, a Wigner crystal with inhomogeneous lat- 
tice spacing, and self-bound clusters of identical fermions. 
These are explained on the basis of numerical few-body 
calculations as well as by deriving simple analytic mod- 
els. A short account of these findings has been published 
recently (27j . 

This paper is organized as follows. Section [XT] intro- 
duces the effective Hamiltonian and discusses aspects of 
the corresponding two-body problem. In Sec. IIII1 we in- 
vestigate the many-body ground-state in the limit of ho- 
mogeneous (repulsive) interactions, i.e., dipoles aligned 
perpendicular to the plane of the ring. Section IIVI fo- 
cuses on the case where the interaction is purely re- 
pulsive but inhomogeneous. Partially attractive interac- 
tions, as occur for small enough angles between dipoles 
and the ring plane, are studied in Sec. [V] Appendix [X] 
contains a concise introduction to the numerically exact 
multi-configurational time-dependent Hartree method 
employed for the computation of few-body ground-state 
properties. 
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Figure 1: (color online) Sketch of the ring-shaped trap (radius 
R, width a± <C R) in the xy plane, with the dipoles moments 
d = d (sin a, 0, cos a) polarized in the xz direction by an 
external field E. 



II. MODEL 

A. Hamiltonian 

We consider a system of N identical particles (bosons 
or fermions) of mass m, confined in a ring-shaped trap- 
ping potential. The ring (radius R) is taken to lie in 
the in/-plane (Fig. [1]). The particles (e.g., atoms with a 
magnetic moment or polar molecules) have a dipole mo- 
ment d = d (sin a, 0, cos a) aligned in the irz-plane by 
an external field, at an angle a to the z-axis. The in- 
teraction between two dipoles in free space is V(r) = 
D 2 (l — 3cos 2 9 r d)/r 3 , where r is the separation between 
the dipoles and 9 r d is the angle between r and d; fur- 
thermore D 2 = <i 2 /47reo for electric dipoles and d 2 /io/47r 
for magnetic ones. Note that for small separations, the 
far-field dipolar interaction is no longer valid, and short- 
range interactions dominate. In order to map out the 
genuine dipolar physics, we assume that short-ranged 
forces can be omitted on the length scale of interest. One 
might imagine tuning them to zero using optical or mag- 
netic Feshbach management; however, it may be impor- 
tant to have short-range repulsion in order to stabilize 
the dipoles. 

We now focus on the limit of a ring potential with tight 
harmonic confinement in the transverse direction. Then, 
for sufficiently weak interaction strength, the transverse 
motion is frozen in the lowest-energy mode, which is 
a Gaussian ip (p) = e~ p / 2aj - / y^a^ of spatial extent 
a± =y/h/mLd± <C R [14]. Averaging out the trans- 
verse degrees of freedom over the reduced density matrix 
/5j_ = \(Po)((p n I > one arrives at an effective ID Hamilto- 
nian [16| 

h 2 N ri 2 

i—1 1 i<j 

where the angle 9i specifies the position of particle i on 
the ring. For R ^> a±, the effective interaction takes 
the form F 1D (0i,0 2 ) = ^cm (6) VW(0), where 6 = {6 X + 
02) /2 is the center-of-mass (CM) angle of the two dipoles 
and i? = 0\ — 02 is the relative angle. In terms of s = 




Figure 2: (color online) Effective interaction: (a) Potential- 
energy surface Vid (01,0a) (here for a — -|). (b) Center-of- 
mass- dependent interaction strength Vcm(O) for a = 0, 0.197T 
and j. (c) Relative-coordinate potential V rc i(i}) [thick line] 
for R/a±_ = 10, including the long-range asymptotics (thin 
line) . 

2R\ sin($/2)|/a_L, the dependence on the relative angle is 
given by 

V re i(0) = v^7r(l + s 2 )e s2/2 erfc(s/\/2) - 2s; (2) 
moreover, 

D 2 

V CM (e) = -^-(1 -3 sin 2 a sin 2 6). 

Note that the effective interaction Vid is inhomo- 
geneous in the angular coordinates: This reflects the 
anisotropy of the underlying 3D interaction, which favors 
configurations where the two dipoles are arranged prefer- 
ably parallel tod, i.e., such that cos 9 r d — sin a sin ( £i+*a ) 
is maximized, with the particles restricted to the ring. 
Thus the CM potential, shown in Fig. [51 has minima at 
= ±-|, which become more pronounced for larger a. 
For 

a > a c = arcsin ~ 0.1967T, 

the potential acquires attractive regions. The potential 
should be viewed as a potential-energy surface which de- 
pends on both particles' coordinates rather than just the 
distance between them, as sketched in Fig. [5] 

The derivation of the effective interaction © assumes 
that the single- mode approximation is valid, i.e., that 
only the lowest transverse mode tpo is occupied. This is 
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equivalent to first-order perturbation theory or, in the 
language of scattering theory, the Born approximation 
(where ipo(p)e lqx replaces the 3D plane- wave states in 
the presence of confinement). This is valid so long as 
the interaction energy is small compared with the trans- 
verse level spacing, ftioj_, i.e., if <C In that 
case, second-order perturbation theory predicts negative 
corrections 0(^ sL ) 2 due to virtual excitation of higher 
transverse levels. For interaction energies of order tiiu±, 
confinement-induced resonances may be expected, simi- 
lar to those for contact interactions under transverse con- 
finement Hi. 



B. Two-body problem: Relative motion 

Let us briefly discuss some properties of the two-body 
relative problem, given by the Hamiltonian 
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Strictly speaking, this is relevant only if CM and relative 
motions decouple - as in the case of dipoles perpendicular 
to the plane of the rin g (a = 0), but also, more generally, 
for a linear geometry [14| . However, it will also provide 
useful insight in cases where the particle spacing is small 
compared with the length scale on which Vcm varies, so 
that one may locally replace D 2 /Aa\ = Vcm(©) at fixed 

e. 

For distances \r\ ^> a±, the free-space dipolar potential 
is restored, so that the Hamiltonian becomes a± indepen- 
dent, 
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Its strength is characterized by the dipolar length 
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At short distances, the transverse average introduces a 
short-range "cutoff" a± . This attenuates the r = diver- 
gence into a peak of order J) 2 /q^.[4^| 

As an illustration, we have diagonalized h TC \ numeri- 
cally and plotted the ground-state wave function i/j re i(?9) 
in Fig. [3] for different rd/R. For bosons, the noninteract- 
ing wave function is the zero- momentum state ip{$) = 
—h=. For rd/R <C 1, the dominant contribution comes 
from the peak of V IC \ at 1? = 0, where the probabil- 
ity amplitude is reduced. With increasing r<2, the mini- 
mum value of ^(O) becomes deeper until it reaches down 
to almost zero (see, e.g., rd/R — 0.4), in which case 
ipi'd) ~ | sin j\ resembles the modulus of the wave func- 
tion of noninteracting fermions. This behavior is reminis- 
cent of bosons with a short-range (contact) interaction. 

For rd/R > 1, the weaker |r| -3 tail becomes effective, 
driving the two dipoles apart until their density is sharply 
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Figure 3: ('color online) Relative wave function ip^) of two 
bosonic dipoles on a ring. For weak repulsive interactions, 
rd/R < 1, the ip resembles that of bosons with contact inter- 
actions. For stronger coupling, rd/R > 1, the wave function 
becomes localized at # = ±7r. For sufficient attraction, the 
bosons form a bound state. Inset: For identical fermions, 
a tightly bound state may form for stronger attraction (here 
r d /R=-4). 



peaked at a distance = ir (i.e, the dipoles are at op- 
posite poles of the ring, |r| = 2R). This is a precursor 
of a crystal-like state. For fermions, the modulus |V>($)I 
would look similar. 

In the spirit of our comment above, let us consider an 
effective attraction, rj, < 0. Figure [3] illustrates that two 
bosons form a bound state, which here becomes local- 
ized for \rd\/R > 0.0 1. This is analogous to the well- 
known bound state of a ID delta interaction. However, 
even identical fermions can form a "p-wave" bound state 
(Fig. O inset), provided the potential is deep enough to 
accommodate a second, anti-symmetric level. Note that 
it owes its existence to the non-zero range ~ a± of the 
potential. 

Let us mention that the relative state for N = 2 
fermions (or any even number) on a ring is a sub- 
tle issue, as can be seen from the Slater determinant 
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The CM state 



has nonzero (angular) momentum ±1, whereas the rel- 
ative orbital, having momentum i, is antiperiodic in 
i9 = 6*i — 62 in order to preserve 27r-periodicity of the 
total wave function. This relates to the fact that any 
even-N Fermi sea is twofold degenerate, and thus will 
either have nonzero momentum or, in the case of a sym- 
metric Fermi-level state e ikF0N + e ~ ikpeN cx cos{k F 6 N ), 
a cosine modulation of the CM wave function, breaking 
the translational symmetry. 



III. HOMOGENEOUS CASE (a = 0) 

Before tackling inhomogeneous interactions, let us first 
focus on the homogeneous case, a = 0. Here the dipoles 
are oriented perpendicular to the plane of the ring; thus, 
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Figure 4: (color online) Pair-distribution function p2(6,0) in 
the homogeneous case a — (4 bosons), for nra — 0.0025, 
0.025, 0.25, 2.5 and 25.5 (with increasing dashing). 
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rotational symmetry is preserved and the interaction is 
independent of 9, V 1D (6 1 ,0 2 ) = D 2 V le i(0i - 2 )/4a\. 
This system therefore has a close analogy to a linear ID 
system, with particle distances Xi — Xj = 2Rsm( 6l ~ 62 ). 
In order to exhibit the essential physics, we first discuss 
ground-state properties for a few particles as obtained 
using the numerically exact multi-configurational time- 
dependent Hartree method (see Appendix EJ) ■ We then 
proceed to discuss analytical models for the limiting cases 
of gas-like and solid-like regimes. 



Figure 5: (color online) Ground-state energy E{nrd) for TV = 
4 dipoles in the homogeneous system (a = 0). Also plotted 
for comparison: the energy of N Lieb-Liniger bosons (with 
g — D 2 /a 2 L ) and the classical limit Ec- 



B. Bosons: Lieb-Liniger-gas regime {nrd <C 1) 

For weak dipolar coupling, r<± <C 1/n, we observe a 
behavior reminiscent of a ID Bose gas with short-range 
interactions. Let us now model this gas-like regime. 

The ID Bose gas is described by the Lieb-Liniger model 

Hi 



A. Ground-state properties 

Since the number density n(9) = (JY S(9 — ftj)) 
Y 'It. is constant for any nrd due to translational invari- 
ance in 9, we focus on two-particle correlations. These 
are described, e.g., by the pair distribution function 
P2(0,0') = £ i7 y<<K0 - 9i)6(ff - 9j)) t which gives the 
probability density of finding one dipole at 9 and a sec- 
ond one at 9'. This is plotted in Fig. 2] for N — 4 bosons. 

For dipole lengths much smaller than the inter-particle 
distance, nrd *C 1, the initially uniform pair distribution 
develops a "correlation hole" at 6 = 6' with increasing 
coupling — indicating that finding two dipoles at the 
same position becomes more and more unlikely — un- 
til the pair distribution strongly resembles that of non- 
interacting fermions (see, e.g., nrd = 0.25) exhibiting 
Friedel oscillations. This gas-like behavior is reminiscent 
of bosons interacting through a short-range or contact 
interaction g5(x — x') (2^ |. an observation similarly made 
in the context of a dipolar gas in a harmonic trap [3uT |. 
Throughout this regime, fermions would remain essen- 
tially noninteracting due to the Pauli principle. 

As the dipolar length becomes comparable with the 
mean inter-particle separation, nrd ^ 1, the long-range 
tail of the interaction becomes crucial. Then the pair dis- 
tribution develops more pronounced oscillations, which 
culminate in well-separated peaks for nrd S> 1. These 
indicate the crystal-like localization of the individual par- 
ticles 
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Xj), 



(3) 



where here Xi = R9i. In the thermodynamic limit 
N, R — > oo (keeping the density n = N/2ttR fixed), the 
properties of the Lieb-Liniger model are completely de- 
termined by the dimensionless parameter [29] 



7 = 



TO 

h 2 n 



>0, 



where a\ = —2h 2 /mg denotes the ID scattering 
length. More specifically, the Lieb-Liniger gas exhibits a 
crossover from a weakly interacting Bose gas for 7 <C 1 to 
a so-called Tonks gas in the "fermionization" limit 7> 1, 
similar to our findings in the gas-like regime. We will 
now delineate a relation between the model parameter 7 
(i.e., g) and the ID dipolar interaction (parametrized by 
rd,a±). 



1. Born approximation for the interaction 

In the spirit of the Born approximation, one 
may identify the coupling constant with the zero- 
momentum Fourier transform of the potential, g w 
V ip (q = 0) = lim^o / dx V 1D {x)e^ x . For V 1D (x) = 



D 2 V Ie i(x/a±)/4a,j_, this yields 



D 2 
2a\ 



ds V, c 



r<2 



2na 2 L 



(4) 



5 



To check this, in Fig. \5\ we have plotted the ground- 
state energy E(nrd) obtained numerically for dipolar in- 
teractions as well as for the Lieb-Liniger model. To com- 
pare the two on an equal footing, we have made use of the 
results for N — 4 bosons based on the solution of the N- 
particle Lieb-Liniger equations (3l| (rather than the inte- 
gral equations valid in the thermodynamic limit), with g 
given by (gj). For 7 < 1, i.e., nr d < 2(na±) 2 ~ 0.01, the 
system is well described by a weakly interacting Bose gas, 
E ~ ^-^-ng. Conversely, for 2(na±) 2 <C nrd <C 1, one 
has 7^-1, and the bosons repel each other so strongly 
that they become virtually impenetrable. For very tight 
confinement, na± — > 0, the bosons are thus practically 
fermionized already at very low coupling, which is the 
case for a purely l/|f| 3 potential [lOj . 

In the strict limit g +00 (a% — > _ ), the Lieb-Liniger 
state \l/ ai is closely connected to that of non-interacting 
identical fermions ^> f via the so-called Bose-Fermi map 

m, 



ai-vO- 



A^ F ; A(xx...x N ) = Y[sgn(xi 



*<3 



Note that for the ground state, ^0 = | "fp | . Consequently, 
only phases differ, whereas all local quantities are iden- 
tical. Thus the pair distribution function (Fig. @| dis- 
plays a correlation hole of size r ~ 1/fcp, kp = irn be- 
ing the Fermi wave number, and the energy approaches 

E -> E W < kF ^ ^ ^e F , where e F = (hk F ) 2 /2m. 
Notice, however, that here the fermionization limit is 
slightly subtle. First, 7 in (|4]) cannot tend to infinity 
so long as nrd *C 1; rather, at 7>l/2(naj_) 2 , nonzero- 
range effects come into play which are no longer de- 
scribed by the Lieb-Liniger model. Moreover, for even 
numbers N — 2,4, . . . , the Bose-Fermi map is not appli- 
cable on a ring [32] , since the Fermi gas has CM momen- 
tum k F = irn for reasons discussed in Sec. Ill Bl Thus, the 
Tonks gas has an energy lower than that of the Fermi gas 
by the CM energy (fifcp) 2 /2Nm. Generally, for any finite 
N, the energy is slightly lower than the thermodynamic- 
limit result E = \Ne F by a factor of (N 2 - 1)/N 2 . 



2. Beyond the Born approximation 

Let us now discuss the validity of the Born approxima- 
tion (0J and how to include effects beyond it. 

In one dimension, the Born approximation holds for 
ID scattering lengths large compared with the scatter- 
ing wavelength, |fcai| 3> 1, and becomes exact in the 
non- interacting limit a\ — > —00. This can be derived 
straightforwardly from the Born expansion of the T ma- 
trix. For low-energy scattering, ka± <C 1, this leads to 
the condition |oi| 3> a±. Since for the dipolar interaction 
we have |ai| ~ a^/r^, this is equivalent to 

r d < a±. 

Thus, the Lieb-Liniger map ((4]) is valid for nrd <C na±, 
tightening the constraint found above. As an illustration, 
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Figure 6: (color online) Effective ID interaction strength g = 
—2h 2 /mai (units of h 2 /maj_) for the dipolar potential Vid 
as a function of raj ax- For comparison: Born approximation 
g « D 2 /a^, and results for a box potential with corresponding 
parameters. 



for N = 4 and R/a± = 10, we have na ± w 0.064, i.e., (g]) 
holds up to 7 ~ l/(naj_). 

So, what corrections for g do we have to expect for 
stronger couplings, > aj_? Although that regime is 
hard to realize experimentally, the question is of theo- 
retical interest. To answer it, we have calculated nu- 
merically the ID scattering length, a\ — — 2h 2 /mg, for 
the dipolar potential Vijj. This is done by solving the 
E = Schrodinger equation for the relative motion, 

+ Vi D {r) tp(r) = 0, and fitting it to the asymp- 
totic low-energy form ip(r) ~ e(|r| — ai) for \r\ 3> aj_. 

The results are shown in Fig. [SJ Clearly, the Born 
approximation remains qualitatively correct up until 
Td/d± ~ 0.5, corresponding to nrd ~ 0.5naj^ ps 0.032 
above. However, as rd/a± — > 0.77..., a resonance ap- 
pears, signifying that g 00, and then g crosses over 
to g — > —00. This resonance explains why, in Fig. [SJ 
the fermionization limit was reached already much earlier 
than naively expected from ([J}. Moreover, it suggests an 
interpretation for energies larger than that of the Tonks 
gas, but still in the gaseous regime: A negative value of 
g indicates the existence of a so-called Super- Tonks state 
|33l | , which exhibits correlations analogous to a hard-core 
potential with nonzero range a\ > 0. Before discussing 
this in more depth, let us illustrate briefly why the reso- 
nance comes about. 

a. Box-potential model To obtain some intuition for 
interpreting the numerical results obtained above, let us 
consider a crude model for the dipolar interaction above: 
a box potential V(x) = v&(£ — \x\) of range £ ~ a± and 
height v ~ D 2 /a 3 L . The ID scattering length is known 
analytically, 



a 1 



1 



coth?7 




Of course, for 77 -C 1, the Born approximation is recov- 
ered, a\jl ~ — 1/77 2 , or g ~ 2£v ~ D 2 /a 2 _. However, at 



7] = T} 



1.2 (i.e., for of order a±), a\ vanishes like 
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ai/£ ~ i(r; 2 — n* 2 ), so that g oc (r<2 — exhibits a 

resonance at the critical point. This is exactly the type 
of behavior encountered in Fig. H2 

To understand this phenomenon a little more deeply, 
it is worth considering what happens to the wave func- 
tion. For r\ <C 1, the (negative) scattering length is much 
larger than the box radius, |ai| ^> £, which reflects that 
the length scale of the wave function is too large to sample 
any details of the potential. In other words, the typical 
energy is so small compared with that associated with 
the potential's structure, h 2 /ma 2 <C h 2 /m£ 2 , that the 
physics is basically shape independent and entirely de- 
termined by the Born parameter. However, when rj ~ 1, 
we have |ai| < I , i.e., the essential change occurs within 
the interior of the box, where the wave function is sup- 
pressed more and more strongly. In the extreme case of 
a hard-sphere potential, n 3> 1, ijj(r) = for all \r\ < £, 
so that the scattering length a\ = I > (in the dipolar 
case, a\ — > +oo). By continuity, a\ must change sign in 
between. 

b. Super-Tonks regime The results above permit 
some analytical statements about the regime where a\ > 
0, i.e., beyond the resonance. There the system can be 
understood as a gas of ID hard spheres (rods) with di- 
ameter 1 = a\. As this resembles a so-called Tonks gas of 
fermionized bosons, with diameter a% — > 0~ , this is some- 
times referred to as the Super-Tonks regime [33l l34j. 

In a homogeneous system, the wave function for na% > 
can be obtained from the fermionized wave function 
(di = 0) via the mapping [35] 

*«i ({%}) - *o- - tt - . 

where ordering Xj < Xj+i — a\ is understood, and the 
total length L entering is replaced by L — Na\. With 
that substitution, the hard-sphere energy can be obtained 
from the fermionized value, E ai ^Q = N(hirn) 2 /6m, 



En 



(1-nai) 2 



(5) 



In principle, one might expect this formula to also 
give a good description of the dipolar ground-state en- 
ergy E(nrd) when inserting a\(rd) as obtained from the 
scattering calculation, at least in the neighborhood of 
na\ = 0. However, in our case this regime is not well 
separated from the long-range one, nrd ^> 1, as explained 
above. Consequently, we do not find Eq. ([5]) to yield a 
reasonable agreement for our parameters. 



C. Localization (nrd 3> 1) 



become classical in the sense that the potential energy 
overwhelms the kinetic energy, so that for nrd — ► oo the 
ground state becomes a position eigenstate \6%, . . . , 0^)± 
(± indicating bosonic/fermionic permutation symmetry), 
where the angles {6j} are determined by minimizing the 
energy 
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This yields a crystal-like, equidistant distribution of par- 
ticles: 9j = + 5, j € {0,...,JV- 1}, which is 
only defined up to an angle 5 by rotational invariance, 



and an energy E y 
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sm 



^| . To find 



the behavior in the thermodynamic limit, one may ap- 

. i—3 

proximate the sum by an integral, sm TT ~ 
2f Jjy^rf^lsin^p 3 ~ 4N 3 /n 3 as TV — > oo. In that 

limit, we thus have Ec/N rj 2nrd x (hn) 2 /2m, an ex- 
pression formally similar to that of an infinite line [10|. 

However, for finite nrd, the particles are not strictly 
localized. Rather, their wave packets are spread out over 
a width w due to their zero-point motion, as is illus- 
trated by the pair distribution function in Fig. [U We 
now estimate w. For simplicity, we consider a linear 
ID system and perform a classical normal-mode analysis: 
Expanding the interaction potential D 2 X)i<j \ x i ~ x j\~~^ 
to second order about the minimum (a lattice with 
nearest- neighbor distance 1 /n) yields effective harmonic- 
oscillator frequencies ujq = -\/l2Z) 2 n 5 /m, corresponding 
to an oscillator width w = yjh/mujc given by 



nw 



12nmD 2 /h 2 



This confirms that for dipole lengths large compared to 
the average inter-particle spacing, nr^ 3> 1, the dipoles 
should indeed be well localized, in agreement with Fig. [4] 
Taking into account the zero-point motion about the 
classical equilibrium positions, we find that the crystal 
energy Ec is corrected by +Nhajc/2 ~ N(hn) 2 /m x 
\fnr d /2. 

In the crystalline limit nrd oo where the behavior is 
nearly classical, the same reasoning applies to fermions. 
For finite coupling, however, due to the Pauli principle, 
the kinetic energy of fermions is inherently higher than 
that of bosons. In this sense, stronger repulsion is needed 
to suppress the zero-point motion compared to bosons, 
so that the crystalline limit is approached more slowly. 



Let us now discuss the limiting case nr^ ^> 1. We 
found that the long-range dipolar repulsion becomes 
dominant here, so that the particles localize. Another 
indicator of this is the ground-state energy of the dipoles 
(Fig. [5]), which deviates from the saturated fermioniza- 
tion regime and changes over to an energy increase pro- 
portional to rd for nrd 3> 1- In this regime, the dipoles 



IV. REPULSIVE INTERACTIONS (0 < a < a c ) 

Let us now consider the case where the tilt angle is 
nonzero such that the effective interaction is still repul- 
sive but inhomogeneous. We first present the ground- 
state properties obtained for the few-particle system, be- 
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fore discussing simple analytic models for the gaseous and 
crystal-like limits. 

A. Ground-state properties 

To bring out the essential effect, we first focus on 
bosons with an inhomogeneity just below the threshold 
value, a = 0.19-7T < a c . In contrast to dipoles aligned per- 
pendicular to the plane of the ring, for a ^ the system 
is no longer rotationally invariant about 9. Thus interac- 
tion effects are reflected in the number density n(8), as 
shown for N — 4 bosons in Fig. [7] For small couplings, 
nrd <C 1, the state resembles that of a Bose gas. The 
density has a slightly stronger weight near the potential 
minima 9 = ±tt/2 for intermediate couplings but tends 
to flatten out for larger values, nrd ~ 1- By contrast, 
for nrd 3> 1, the density profile develops N = 4 distinct 
peaks, revealing a crystal-like localization of the parti- 
cles. Note that the peak positions are not equidistant, 
but rather slightly displaced toward 9 = ±7r/2 due to 
the pronounced minima of Vcm-[50] 

The transition from a gas-like to a localized state is also 
clearly visible in the pair distribution function p2{9,9') 
shown in Fig. [7) For nrd "C 1, two-body correlations are 
absent in the mean-field regime, pi(9\, #2) ~ p{d\)p(92), 
or are limited to Friedel-like oscillations around the cor- 
relation hole near 9\ =92, reminiscent of a fermionized 
Bose gas. For nrd ^> 1, the pinning of the particles to 
individual peaks is clearly discernible. 

A similar trend would be observed for fermions: For 
nrd "C 1, we find a gas-like state, which however is only 
weakly interacting owing to the Pauli exclusion principle. 
With increasing coupling, nrd > 1, the state crosses over 
into a crystal-like one characterized by localization of the 
particles in the classical minimum-energy configuration, 
just as in the Bose case. Due to the Fermi energy, that 
crossover is much smoother than for the bosons. 



B. Inhomogeneous Lieb-Liniger regime (nrd <C 1) 

Let us dwell on the parameter range where the ground 
state is gas-like. In the homogeneous limit a — ¥ 0, we 
established that here the system is well described by the 
Lieb-Liniger Hamiltonian ([3]), with the ID interaction 
strength g = D 2 /a 2 l (for <C a±). By extension, for 
a / 0, the system should resemble a Bose gas with an 
inhomogeneous interaction g ( Xl ~^ X2 ) 8(x\ — X2), where 
g(x) = D 2 (l — 3 sin 2 a sin 2 ^)/a 2 _ is the Fourier trans- 
form of V\d at zero relative momentum. 

Although this inhomogeneous Lieb-Liniger Hamilto- 
nian system is not integrable, some insight can be 
gained by assuming slow variation of g(x) compared 
with the local correlation length. In this case, a local- 
density approximation [36] may be applied, which con- 
sists in assuming the equation of state of the homoge- 
neous Lieb-Liniger system, /i = f g (n), to hold locally, 





n r d = 0.0025 






0.18 






1.28 — 






5.1 






\ A 25.5 


.... f:, ; A 





















-1 -0.5 0.5 1 

e/71 




Figure 7: (color online) Inhomogeneous case (a = 0.197T, N = 
4 bosons): Density profile n(9) [top] and pair distribution 
function p 2 {6i,6 2 ) [below] for nr d = 0.025, 1.28, 25.5 (left to 
right). 



p = fg( x ) [n(x)]. Here p, denotes the chemical potential 
as a function of number density n[x); moreover, for the 
Lieb-Liniger system, f g (n) can be obtained directly from 
the ground-state energy, which is known in terms of the 

dimensionless function e(i) = E/N ( -^- (HH^]. Note 
that the chemical potential is x independent, as local 
equilibrium is assumed to hold. Carrying out the local- 
density approximation numerically for the full range of 
7 is beyond the present scope. However, considering the 
borderline case also provides some insight into the impact 
of the inhomogeneity. 



1. Mean-field regime, y(x) <C 1 

For locally weak interactions, 7 <C 1, the homogeneous 
equation of state reads // = ng, which yields a profile 

n ( x ) = = "(0) 
{ ' g(x) 1 - 3 sin 2 a sin 2 |' 

Remarkably, this is independent of 7, and thus of nrd, 
since pocg, that is, kinetic-energy contributions are neg- 
ligible in the homogeneous equation of state. This holds 
only for large enough particle numbers N 3> 1 and not 
too large variations of g(x); otherwise the kinetic pres- 
sure will smear out the density gradients on the scale of 
the coherence length or smaller. 

It is worth mentioning that the local-density approx- 
imation for 7 -c 1 corresponds to the Thomas-Fermi 
approximation, which neglects the kinetic energy in the 
Gross-Pitaevskii equation for the mean-field orbital <fi(x) 
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[H, Ch. 6], 



2m 



4>"{x) = [m - n(x)g(x)] <p(x) 



0. 



While the strict mean-field limit assumes all bosons 
to be condensed into a sing le orbital, * = <f)® N /N N / 2 , 
even for weak repulsion a deviation from a fully uncor- 
rected state occurs. Conveniently, this can be quan- 
tified, e.g., in the second-order correlation function, 
g2(x,x) = p2(x, x)/n(x) 2 . In the homogeneous case, 
52 = 1 - + 0( 7 ) [H. Replacing 7 = mg/nh 2 

by 7(2;) — mg(x,x)/n(x)h 2 ~ mg 2 (x, x)/h 2 [i yields 
.92(2^, x) ~ 1 — 7 V7(0)(l — 3sin 2 a sin 2 -^), amounting to 
a pronounced suppression of the pair distribution pi(9, 9) 
near 9 = x/R = 0, as observed in Fig. [Jj 



2. Tonks-gas regime, y(x) ^> 1 

For stronger coupling, the density modulation (Fig. [7J 
nrd = 1.3) is attenuated until, in the fermionization limit 
7(2;) ^ 1, n(x) approaches a constant. This is clear since, 
by the Pauli principle, identical fermions do not sense 
any contact interaction and hence no spatial modulation 
of g{x). Formally, this can be seen from the chemical 
potential in the Tonks regime, p, = (fom) 2 /2m, which 
becomes independent of g. Since fi is spatially constant 
by assumption, so is n(x) = n(0). 

For large but finite g, pair correlations are actu- 
ally not fully suppressed, 32 — + |(7r/7) 2 . Since 
7(2:) ~ mg(x)/nh 2 , this implies fairly strong residual 
pair correlations near the potential minima 9 = ±7r/2, 

52(2:) ~ 37(0)^ ~ 3 sin 2 a sin 2 j^J . Moreover, as dis- 
cussed for the homogeneous case, the limit 7(2;) ^> 1 
corresponds to fermionization only if nrd <C 1, locally. 
Instead, for increasing nrd, we find that the incipient 



fermionization gives way to a state where the nonzero 
range of the repulsive potential becomes relevant, until, 
eventually, a crystal-like localization occurs. 



C. Inhomogeneous crystal-like regime (nrd S> 1) 

We now present analytical models for the crystal-like 
limit nrd 3> 1. Similarly to the case a = , here each 
dipole is pinned to a position 9i determined by the ex- 
trema of the classical energy 



E< 



D 2 v 1 

i<3 



o • 2 -2 s t +e 
6 sin a sin — 15- 



sm 



In this case, though, the equilibrium configuration is not 
simply an equidistant lattice. Rather, owing to the com- 
petition of potential energy of the relative motion (fa- 
voring maximum distance) and the CM (alignment near 
9 = ±7r/2), the dipoles are localized at discrete positions 



which are distributed inhomogeneously on the ring. In 
principle, the distribution of {9i} can be found by nu- 
merical minimization of Eq- However, to obtain some 
analytic insight, we will do this for one of the simplest 
nontrivial cases, N = 4. Complementing this few-body 
viewpoint, we then extend the local-density approxima- 
tion above to the strong-coupling regime. 



1. Few-particle model 

Let us consider the special case N = 4. By symmetry, 
we can write §i = ±(? ± </>), i = 1,...,4. The one- 
dimensional potential-energy curve 



Ec(4>) 
D 2 /AR 3 



1 — 3 sin a 



sm 



1 



+ (1 



3 sin 2 a sin 2 1 



consists of three contributions: 



1. The first term represents the weak repulsion of a 



pair at, say, 

(e = tt/2). 



• ± (j), which has minimum CM energy 



2. The second term represents the energy of a con- 
figuration of type ±(^ + </>), which has maximum 
repulsion (O = 0). 

3. The remaining term is that of two diametrically 
opposed dipoles, e.g., at ±^ + 0; here only the 
slowly varying CM potential contributes. 

The total potential consists mainly of (1.) a steep poten- 
tial barrier at (j) — > (repulsion of a pair near ±7r/2) and 
(2.) a confinement for larger (f> due to inter-pair repul- 
sion, adding up to a potential well for <fi. Neglecting the 
third term, the potential minimum can be found to be at 



0* (a) — arctan yl — 3 



sin 2 a. 



(G) 



Near a = 0, this reproduces the uniform lattice spac- 
ing 2n/N, with a very smooth correction for a > 0, 



2</>* (a) 



3„,2 



a . It is only very close to the critical 



angle that a noticeable inhomogeneity appears, 0* (a) ~ 

V^l — 3 sin 2 a ~ 8 5 / 8 $a c — a. At a — > a~ , this goes 
to zero very abruptly, indicating an attempted clustering 
as the minimum energy vanishes. Including the slightly 
deconfining third term above, the minimum separation 
angle is marginally shifted toward higher (p. Moreover, 
although we chose N = 4 as a special soluble case, the ba- 
sic competition between bunching toward 9 = ±7r/2 and 
keeping a minimum inter-particle distance will persist for 
the case of more than a single pair, N > 4. 

Let us comment on the limit a — > a c . The non- 
analyticity of <fi*(a) at a c hints at the fact that, at least 
for nrd ~> 00, there is a phase transition when a crosses 
the critical angle. In fact, we will see in Sec. |V]that for 
a > a c bound states appear which are tightly localized 
at 9 = ±7r/2 for nrd oo- For a < a c . however, the 
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two dipoles at, say, tt/2 cannot get arbitrarily close due 
to the quantum-mechanical kinetic energy neglected so 
far, contrary to what we would expect based on 



2. Local- density approximation 

Complementary to the few-particle model discussed 
above, let us consider the limit N — > oo. If the CM po- 
tential varies on a length scale much larger than the aver- 
age inter-particle distance, we may assume that, locally, 
the system behaves similarly to a homogeneous ID sys- 
tem. Then the local-density approximation may be ap- 
plied, i.e., that the equation of state of the homogeneous 
system, [i = f Td (n), holds locally fi = / rd(e) [n(0)]. 
Here n(0) is the average number density around 0, and 
r d (Q) = r d x (1-3 sin 2 a sin 2 ©). 

If we take the chemical potential of the linear homo- 
geneous system for nr d — > oo, fi — 4£(3)7i 2 n 3 r ( j/m ficfr . 
we find 



n(6) oc 



1 



V 7 ! - 3 sin 2 a sin 2 6 



This density variation is rather smooth due to the strong 
dipolar repulsion. For smaller a, the density is practi- 
cally homogeneous, with tiny humps near O = ±7r/2. For 
a a~, though, the density becomes markedly peaked 
and tends toward (5(0 =p -|), indicating an instability to- 
ward clustering beyond a c . Of course, the local-density 
approximation is valid only for density variations slow 
compared with the local coherence length: Close to the 
critical angle, n(0) is smeared out on the length scale 
given by the zero-point motion. 



V. PARTIALLY ATTRACTIVE INTERACTIONS 

(q c < a < « C 2) 

For tilt angles a > a c , the interaction potential is no 
longer purely repulsive but acquires attractive regions for 
CM angles |0| = § ± 4f (Fig. 0), where 



A0 = 2 arccos 



Ax 
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Figure 8: (color online) Clustering of N = 4 bosons (a = 
0.2n): Density profile n(6) [top] and pair distribution function 
P2{0 1 ,9 2 ) [below] for nr d = 0.025, 0.802 and 0.807 (left to 
right). 
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Figure 9: (color online) Clustering of iV = 3 fermions (a = 
0.22-7T): Density profile n{9) [top] and pair distribution func- 
tion p2(6>i,6> 2 ) [below] for nr d = 0.57, 2.14 and 2.20 (left to 
right) . 



Note that for a — > a+ , this attractive segment is vanish- 
ingly narrow, A0 ~ 2 7 / 4 ^a — a c , whereas for a — > \, 
it makes up almost two thirds of the total circumfer- 
ence, A0/7T — !• 0.608 .... With increasing coupling, the 
competition between repulsion and attraction leads to an 
intriguing crossover from gas-like states toward clustered 
"droplets", which we will first discuss on the basis of few- 
body simulations, before giving an analytic description 
of the cluster states. 



A. Crossover from gas-like to clustered states 

Let us first illustrate the crossover from weak to strong 
coupling by means of the density n(8), shown in Fig. [5] 
(N = 4 bosons, a = 0.2tt) and Fig. M (N = 3 fermions, 
a = 0.227r). In both cases, the densities go over from a 
homogeneous profile, characteristic of a gas-like state, to 
one with two sharp peaks, indicating cluster formation. 
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Figure 10: (color online) Ground-state energy E(nrd) in the 
case of partial attraction: N = 3 fermions ( — , a — 0.227r); 

N = 4 bosons ( , a — 0.2ir, 0.227T, and n/2 from top to 

bottom). 



It is noteworthy that in these examples, the transi- 
tion to the localized clustered state is very sharp (e.g., 
at nrd R3 0.8 for the bosons). This is also reflected in the 
ground-state energy E(nrd), shown in Fig. \W[ In fact, 
whereas the energy eventually falls and becomes nega- 
tive at high enough nrd for any a > a c , the behavior 
for small coupling crucially depends on a: Closer to the 
critical angle, the bosonic energy increases for small cou- 
pling, and it is only beyond a critical value of nrd that it 
starts decreasing rather abruptly. By contrast, for larger 
a, E(nrd) is a monotonically decreasing function. To 
understand that competition between repulsion and at- 
traction, note that — for m\i — > — the (bosonic) energy 
behaves as 



E 



N(N — 1) 



<F CM (e)>0 <FrelW> 



(7) 



N 



"5(0) (1 



sin a 



where the average is over the (uniform) noninteracting 
CM and relative states. Since g(0) > 0, the energy in- 
crease for small nrd is positive for 



a < &c2 = arcsm 



0.304tt. 



In other words, for tilt angles less than a c , the potential is 
predominantly repulsive. Initially, the delocalized bosons 
will thus undergo an analogous evolution as in Sec. IIV1 
e.g., behave like an inhomogeneous Bose gas, which is 
predominantly repulsive (a < a C 2) or has a net attraction 
(a > a C 2). For identical fermions, the interaction-energy 
shift per particle will be negligible compared with the 
Fermi energy for all a. 

This makes clear that the transition to a clustered 
state also strongly depends on a. Beyond a C 2, there is 
a smooth crossover from an attractive gaseous state to 
a clustered one (Fig. [TU]) at rather small coupling. By 
contrast, for only slightly overcrtical a £ (a c ,a C 2) - e.g., 
a = 0.27T in Fig. [5]- the transition takes place at much 



larger coupling and is more subtle. With increasing nrd, 
the repulsive gas-like state first goes over into a transi- 
tional state with strong fragmentation (cf. nrd — 0.802), 
which can be thought of roughly as a Mott-insulator- 
type state | y> t) with N/2 dipoles at each of the poles, 
6 = ±7r/2. This can be best seen from the pair distribu- 
tion p2(9, Q'), Fig. Reenter), which upon closer inspection 
confirms that a measurement of one particle at, say, 7r/2 
yields a measurement of (-^ — 1) particles at tt/2 and ^ 
at — 7r/2. There is a sharp avoided crossing to a clus- 
tered state {nrd = 0.807), with all N particles localized 
at either ±7r/2. By parity symmetry, the exact ground 
state is actually a superposition or Schrodinger-cat state, 
^= (I-/V, 0) + |0, N)), as is evident from the pair distribu- 
tion function in Fig. [5] A small symmetry-breaking per- 
turbation will likely cause a collapse of the superposition. 

For fermions, the critical interaction strength needed 
for cluster formation is higher (Fig. ITU]) . This is because 
of the Pauli principle: The exchange hole discernible in 
P2 (Fig. strongly reduces the average attraction be- 
tween two fermions. Moreover, due to the Pauli pressure, 
the fermions are more spread out spatially and require a 
stronger attraction to be squeezed into the narrow zone 
of attraction. 

To obtain a deeper understanding of cluster formation, 
including its dependence on the particle number and on 
a, we will now develop a simplified analytical model that 
can reproduce its basic features. 



B. Simple cluster model: Homogeneous case 

As a first step, let us consider the case of a homoge- 
neous system with purely attractive interaction. This is 
interesting in its own right, since clustered states are not 
limited to a ring geometry. Moreover, as we will see be- 
low, this also allows for a qualitative understanding of 
the inhomogeneous case. 



1. Bosons 

Consider N bosonic dipoles interacting with a homo- 
geneous attractive contact interaction g < 0. (On the 
ring, we may identify this with the value of g(x) in the 
minimum x /R = ±7r/2, g = D 2 (l - 3 sin 2 a)fa]_ < 0.) 
Let us assume all bosons to be localized at some xn, 
spread out over a length L, e.g., by making a varia- 
tional ansatz with all particles in a single orbital, 0l(x) = 
e -(x-x ) 2 /2L 2 j _ xheil; p rov ided L > a±, the inter- 
action energy scales as E[ nt ~ N 2 g/L, whereas E^ ~ 
NH 2 /mL 2 . For g < 0, the total energy E(L) has a local 
minimum, which yields the size of an equilibrium cluster, 



N' 



a i 



2f? 
mg 



(8) 



The characteristic length in the two-body case is, of 
course, the local ID scattering length eti > 0, and larger 



11 



N lead to a shrinking of the cluster. The cluster's binding 
energy is then 



E N A 



-N' 



mat 



ma . | 



which diverges quadratically with the coupling strength. 

Note that in the language of the Lieb-Liniger model, 
this can be understood as the well-known "bright-soliton" 
ground state of attractive ID bosons (iol l4l|: 



E = - 



N(N 2 -1) h 2 



6 



mat 



Since E/N oc N 2 , the binding energy is unbounded in 
the thermodynamic limit N = hL — > oo, and the state is 
stable only for a small enough number of particles. 

In free space, this cluster size would be extremely large 
for weak attraction. On the ring, however, for the cluster 
to form, it must fit into the attractive regions of length 
Ax = 2i?arccos(l/-\/3sina) around 6 — ±7r/2, 



L < Ax. 



(9) 



This yields a rough estimate for the critical interaction 
strength necessary for cluster formation, 
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■C-rMi 

\ v 3 sin a J 



1—3 sin a I 



(10) 



which qualitatively reproduces the parameters necessary 
for cluster formation found numerically (27l |. 



2. Fermions 

We now discuss the fermionic case. As above, we as- 
sume the fermions to be in an essentially noninteracting 
state, but with a localization length L. Then the kinetic 
energy, E kin = N(hk F ) 2 /6m = N 3 n 2 h 2 /6mL 2 , with 
kp = nN/L, has the same structure as in the bosonic 
case except for a factor of order iV 2 due to the exclusion 
principle. By contrast, the interaction energy has a more 
subtle form since the nonzero range of the potential is 
crucial. Setting again D 2 (l — 3sin 2 a)/a 2 _ = g < and 



p 2 (x,x') = (f ) g 2 (x - x') 



we can write 



E\ n t = 



N 2 g 



^ L 92(r)V Iel (— ) . (11) 
o *a,± \ a ±. 



For an ideal Fermi gas, 
92 (r) = 1 



sin(^r) 



7Vsin(£ r) 



JV»1 



sin(/cirr) 



There is no general analytic formula for the integral (fTTjl , 
but we can make statements about limiting cases. First, 



as L <C a± tends to zero, the integral gives LV Ie i(0)/2a±, 
so 



N*g 
2o_l 



Being bounded, the interaction energy is thus completely 
outweighed by the kinetic energy (oc L~ 2 ), which pre- 
vents the collapse toward high-density clusters. On the 
other hand, for large sizes, L Na±, the dominant con- 
tribution comes from the intermediate zone a± <C r <C 
L/ttN where V Te i(s) ~ 4/s 3 and gz(r) ~ \{kpr) 2 , and 



N3 _ al ln(L/nNa ± ) [L 



This falls off slightly faster with L than the kinetic en- 
ergy. Hence, fermions cannot form a bound state with 
arbitrarily large size. This reflects the fact that, for low 
densities, the attractive interaction energy will be largely 
canceled by the exchange term. However, binding is pos- 
sible for an intermediate regime, a± <^ L <C Naj_. Here 
the density is high enough so that the direct (Hartree) 
interaction term is large, ~ 1, but not so high for 

binding to be suppressed by the kinetic pressure: 



E u 



N 2 
2L 



—g (a± < L < Na±) 



It is only for such moderate sizes L that clusters may 
exist. Under this assumption, the minimum energy is 
found straightforwardly for 



Lps — Nax. 
3 



(12) 



By the Hartree approximation, this has the same struc- 
ture as for the Bose cluster, but here the cluster becomes 
larger for higher particle number N because of the Pauli 
principle. Note that this is valid only under two con- 
straints: 

(i) L <C Naj_ for the Hartree term to be dominant, i.e., 
a\ <C aj_ or |r ( ;( : |)| 3> a±. Minimizing E(L) numerically, 
we find the onset of clustering at > a±, which 
relaxes the condition. 

(ii) |^| <^ h 2 jma\ for the strictly ID description to 
hold. This is equivalent to a\ 3> aj_, the complementary 
condition to the Hartree limit (i). However, since the 
latter criterion can be relaxed considerably, it is plausible 
that Fermi clustering may be found in a regime where the 
physics is no longer strictly ID. This may require further 
study, going beyond the bare single-mode ID description. 

The discussion so far assumed a ID system with a ho- 
mogeneous attraction. On the ring, the cluster must fur- 
thermore fit into the attractive regions, L < Ax ©• This 
yields the estimate 



nr d 
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arccos 



|1 — 3 sin a I 



(13) 
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which essentially differs from the bosonic one (JTUJ) by 
a factor of order iV 2 . This turns out to reasonably re- 
produce the curve nr<i{pt) obtained from our numerical 
results (not shown here). 



C. Inhomogeneous case 

Let us now develop a model that goes beyond the in- 
tuitive criterion @ by including the CM dependence of 
the energy. 



1. Bosons 

So far, in assuming a net attraction g < 0, we have 
limited ourselves to a cluster state already tightly local- 
ized in the potential minimum of V"cm- More generally, 
a CM distribution smeared out about — ±^ will shift 
the average interaction energy upward, possibly even to 
positive values, i.e., E- m t ~ ^(g(X = RQ)) with 

(g(X)) = g(R^) + 3 (0)3 sin 2 a(sin 2 (6 -l))=g + Sg. 

For a droplet much smaller than the ring radius, L <C 

R, we have (sin 2 (9 — ■§)) — | (4) ■ This regime ap- 
plies to a tightly bound cluster with a > a c or nrd well 
above the critical value. Then the total energy has an 
interaction term oc L x N 2 g(0) sin 2 a/R 2 corresponding 
to a linear confinement of the CM, in addition to the 
attractive term known from the homogeneous case. For 
convenience, the total energy can be cast in the form 
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happens to give the leading order of the attractive-zone 
width Aa; for a — ¥ a c , and differs only marginally even 
for larger a. For d%/N <C Sx, that is, a broad attractive 
zone which would easily accommodate a cluster with free- 
space binding length d\/N, the solution is 



ai/N 




Sx 



To leading order, this recovers a droplet with the free- 
space binding length, but slightly compressed due to the 
confining term, repelling the CM from the walls of Vcm 
as the CM wavepacket spreads out beyond the potential 
minimum. Conversely, in the case where the free-space 
cluster would have a size strongly exceeding the attrac- 
tive region, d\/N 3> Sx, the droplet size is indeed much 
smaller than that: 



L~ ^/Sx 2 d 1 /N. 



In this limit, the total energy is clearly positive, and the 
object is held together by the inhomogeneity — the CM 
confinement — rather than the attractive mechanism re- 
sponsible in the homogeneous case. 

Note that our above estimate required L <C R, so this 
puts a bound on how weakly confined the dipoles may 
be. For larger droplet size L > R, the CM-confinement 
term Sg will become weaker until eventually, as L 3> R, 
it will saturate to (sin 2 (0 — ■§))—> 5. In that limit, 
the minimum value g will simply be replaced by the uni- 
form average (g) = g + ,g(0)| sin 2 a = (1 — | sin 2 a)g(0). 
This is exactly equivalent to the result obtained the ho- 
mogeneous case ([?]): A bound state can only exist for 
(g) < 0, i.e., in the parameter regime of average attrac- 
tion, a > a C 2- Its binding length is then given by the 
analog of (JSJ, with d\ H> —2h 2 /m(g) significantly larger 
due to the weaker binding. Since L > R by assumption, 
this can describe a cluster state with L < Ax at the 
utmost for large a w n/2. 



2. Fermions 

For fermions, assuming a local-density approximation, 
the inhomogeneity will have a similar effect, namely, to 
smear out the position of the cluster about the mini- 
mum, so that the average direct-interaction term is in- 
creased, (g(X)) = g+Sg. Under the assumptions made in 
Sec. IV Bl this leads to the same results as in the bosonic 
case but with d\/N replaced by Nd\\ For Nd\ -C Sx, 
the cluster size L is slightly compressed compared with 
the homogeneous estimate, Ndi. In the opposite limit, 
L ~ \/Sx 2 Ndi is somewhat larger than the attractive 
zone, corresponding to a repulsively bound object. Of 
course, these estimates are valid only for moderate sizes, 
aj_ <C L < Na±, as discussed above. 



D. Two- mode description 

After having derived a simple model for the single- 
cluster state, let us now take a broader perspective on 
the different phases we found in the regime of strongly in- 
homogeneous, partially attractive interactions. For sim- 
plicity, let us focus on bosons; the fermionic case can be 
discussed analogously. 

• For small enough coupling, the state is always that 
of a weakly interacting gas. For bosons, this is 
a Bose gas with net attraction for a > a C 2 (net 
repulsion otherwise), distributed coherently over 
the regions x± = ±irR/2. If the single-particle 
orbital occupied by all bosons can be written in 
terms of two modes localized at x±, fo(x) = 
^■7= [ip+(x) + <p-(x)], each created by the field op- 
erator a±, then the many-body ground state takes 
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the form 



N 



V2 



|0). 



(14) 



• For sufficiently large nr^ , a cluster state is formed 
for any a > a c . This may be written as a cat state 



1 



\0,N)) = 



1 



(a»+a»y\0). (15) 



• Furthermore, for small enough a > a c , there is a 
transitional regime for intermediate coupling where 
the ground state has a Mott-insulator-like charac- 
ter: That is. half of the dipoles are localized near 
9 = +7r/2 and the other half near —ir/2 (taking 
N e 2N for simplicity), and the wave function ap- 
proximately has the structure 



I 2 : 



1 

2 ' 



,t 



N/2 



N/2 



(16) 



Expanding the many-body Hamiltonian H in terms of 
the two-mode states |JV+,iV_), one obtains an effective 
Hubbard-like Hamiltonian 



Hj,u = - J ^ a\a s > + " s (" s - 1 ) 

(s,s') s 



(17) 



where J = ( L P-\ 2 l^\ ( P+) denotes an effective tunnel cou- 
pling, U = (ip s , ip s \V(x, x')\<p s , ip s ) is the on-site interac- 
tion, and n s = aja s (s = ±). As in the derivation of the 
Hubbard model in the context of external (periodic or 
double- well) potentials [42| , this implicitly assumes tight 
binding, which requires both tunnel coupling and on-site 
interaction to be negligible compared with the excitation 
energy of higher "bands" (i.e., excited states pertaining 
to the doublet <p±). Moreover, we have suppressed an 
overall kinetic-energy offset ~ Nh 2 /mL 2 and an off-site 
interaction U±n+h-, falling off as (a±/R) 3 . 

It is clear from the considerations in Sec. IV Bl that the 
orbitals ip± (and thus J, U) are not known a priori but 
rather depend on N and D 2 . In this sense, diagonalizing 
the effective Hamiltonian Hj t u only yields self-consistent 
information on the ground state. Nonetheless, mapping 
the phases of the two-site Bose-Hubbard model [43| to 
our system may provide some useful insight: 

• For weak on-site interaction, N\U\ -C J, the ground 
state is a delocalized superfluid state (IT41 . maxi- 
mizing the coherence (er, <x_). 

• In the limit of strong attraction, \U\ J, the 
ground state localizes on a single site, thus maxi- 
mizing the absolute interaction energy ~ \U\N 2 /2. 
By symmetry, the state will be a superposition of 
all single-site states, corresponding to the cluster 
state ([I5j). 



• For strong repulsion, NU ^> J, the energy is dom- 
inated by the on-site repulsion. That is minimized 
by a fragmented or Mott-insulator state ([TO]) with 
interaction energy ~ U(N/2) 2 . 

This indicates that the transition from a fragmented to 
a clustered state should take place when U ~ (g(X))/L 
switches from positive to negative. Clearly, this cannot 
occur for a > a C 2 , in which case there is a direct crossover 
from gas-like to clustered behavior. Near the critical an- 
gle in turn, a c < a < a C 2, U changes from positive to 
negative at a±/N < Ax, i.e., for nr^ given by the crite- 
rion (fT0| . Below that critical coupling, the "repulsively 
bound" cluster state (see Sec. IV C 1[) is an excitation on 
top of the insulator- type state, whereas for higher nr^, 
the attractively bound cluster is the ground state. 

VI. CONCLUSION AND OUTLOOK 

In this article, we have shown that dipolar particles 
confined to a quasi- ID ring exhibit inhomogeneous gas-, 
solid- and (clustered) droplet-like ground states, due to 
the underlying anisotropy of the 3D dipole interaction. 

The solid-like behavior is essentially independent of the 
transverse confinement, but differs from that in a lin- 
ear geometry by a non-equidistant lattice-site distribu- 
tion whenever the dipoles are inclined toward the plane 
of the ring (a ^ 0). By contrast, especially the Bose- gas- 
like states depend crucially on the fact that the trans- 
verse confinement length cuts off the |r| -3 divergence of 
the dipole interaction. This leads to a whole class of 
states with a tunable interaction parameter, as described 
by the (inhomogeneous) Lieb-Liniger model, rather than 
merely a hard-core gas. By the same token, the appear- 
ance of attractive interaction regions for dipoles oriented 
sufficiently close to the plane of the ring leads to bound 
states clustered around the "poles" of the ring, rather 
than collapse. 

Our numerical results are presented for few particles; 
however, complementary analytical models support an 
extrapolation to many-body systems. All results are 
based on a Born approximation for the effective ID in- 
teraction, i.e., on the occupation of a single transverse 
mode. An interesting subject of future investigation may 
thus be how an interaction-induced coupling of several 
transverse modes affects these results. Moreover, our cal- 
culations provide motivation for experimental studies as 
well as theoretical extensions to dipolar systems in more 
general curv ed g eometries (including, e.g., coupled quasi- 
1D systems |44l45j). 
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Appendix A: Computational method 

We use the numerically exact multi-configurational 
time-dependent Hartree method [46], a quantum- 
dynamics tool which has been applied to few-body sys- 
tems of identical bosons as well as mixtures (see, e.g., 
[47} for details). Its principal idea is to solve the time- 
dependent Schrodinger equation iMf(t) = H^(t) as an 
initial- value problem by expanding the solution in terms 
of direct (or Hartree) products $j = ipj 1 ® • • ■ S3 fj N ■ 

9(e ] t) = Y t A i {t)^(0;t). (Al) 
i 

Both coefficients Aj (t) and the single-particles basis func- 
tions ifj(t) are determined from the Dirac-Frenkel vari- 
ational principle (6V\[ihd t - H(t)]$(t)) = [46]. This 
leads to a coupled system of Schrodinger- and mean-field- 
type equations for Aj and (pj, respectively, which is in- 



tegrated numerically upon discretization. Note that the 
expansion coefficients inherit the (bosonic or fermionic) 
permutation symmetry of the wave function, Ap(j\ = 
(±l)inv(P)^._ Moreover, in our case of particles on a ring, 
the single-particle wave functions ipj € L 2 (— 7r, it) are pe- 
riodic under translations 9 H> 9 + 2tt, which is ensured 
by expanding them in an appropriate basis set. 

Although designed for time-dependent simulations, it 
is also possible to apply this approach to stationary 
states. This is done via relaxation, i.e., by applying the 
non- unitary imaginary-time evolution operator e~ Hr . As 
t — > 00, this exponentially damps out all contributions 
as e -( E m-E a )r/h exce pt f rorn the ground state. In 
practice, one relies on a more robust scheme termed im- 
proved relaxation (48| . Here (*B\H \ty is minimized with 
respect to both the coefficients Aj and the orbitals ipj. 
The effective eigenvalue problems thus obtained are then 
solved iteratively by first solving for Aj with fixed or- 
bitals and then 'optimizing' ipj by propagating them in 
imaginary time over a short period. That cycle is then 
repeated. 
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